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We generalize the Moyal equation, which describes the dynamics of quantum observables in phase 
space, to quantum systems coupled to a reservoir. It is shown that phase space observables become 
functionals of fluctuating noise forces introduced by the coupling to the reservoir. For Markovian 
reservoirs, the Moyal equation turns into a functional differential equation in which the reservoir’s 
effect can be described by a single parameter. 
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I. INTRODUCTION 

In the study of quantum systems, researchers employ 
different methods to achieve their goals. For instance, 
the dynamics of a system can be described using the 
Schrodinger equation for a state or the Heisenberg equa¬ 
tion for operators. Another example is to use the Wigner 
function W{q,p) to describe the state by a function in 
phase space rather than by a wavefunction ip{q) or a den¬ 
sity matrix p. 

The Wigner function [1, 2], which for a one¬ 
dimensional system takes the form 

w ( q ,p) =-^ s [p]( q ,p) (1) 

s\p](q,p) = Jdq' <q + W\p\q-W) e~ Wp/h , (2) 

has become a very popular tool to visualize the state of 
a quantum system and to compare it to classical sys¬ 
tems. To evaluate mean values of observables A in phase 
space one needs to introduce Weyl symbols S[A](q,p), 
which represent operators on Hilbert space by functions 
on phase space [3, 4], To keep the notation concise we 
will write A(q,p) instead 5'[A](g,p), i.e., the symbol has 
the same notation as the operator but without a hat. 

The dynamics of the Wigner function is described by 
dtW = {H,W}m, where H is the Weyl symbol of the 
Hamiltonian. The Moyal bracket of two Weyl symbols 
A, B is defined as [5-8] 

{A,B} m = ^-(A*B-B*A). (3) 

in 

The star product * between two symbols is a representa¬ 
tion of the non-commutative product between operators 
in phase space [9-13]. For our purpose it can be written 
as 

A(q,p)* B(q,p) = A(C)B(q,p) = B(TZ)A(q,p) (4) 
A(C) = A [q + i^d p , p-i^d^j (5) 

B(Jl) = B (^q - i^d p , p + . (6) 

In the limit h — > 0, the Moyal bracket turns into the 
classical Poisson bracket and the dynamics of the Wigner 
function is described by the Liouville equation [5]. 


The dynamical equation for the Wigner function is 
the analogue of the Schrodinger equation in phase space. 
Likewise, the analogue of the Heisenberg equation of mo¬ 
tion for operators A is the Moyal equation 

dtA = {A, H}m- (7) 

By combining the two possibilities (Hilbert space and 
phase space) to represent a quantum system with the 
two possibilities to describe its dynamics (Schrodinger 
picture and Heisenberg picture), we thus arrive at four 
different ways to study its evolution. 

The dynamical equations above apply to perfectly iso¬ 
lated quantum systems. However, for a realistic descrip¬ 
tion of experiments, the coupling to the environment has 
to be taken into account. This is accomplished by consid¬ 
ering open quantum systems [14, 15], where the system of 
interest is coupled to another system (the reservoir) that 
cannot be accessed in an experiment. Loss of information 
about correlations with the reservoir will introduce deco¬ 
herence and noise to the system. In a similar way, the 
measurement process can be generalized by coupling the 
system to a second system that represents the detector 
[16]- 

In open quantum systems, the four dynamical equa¬ 
tions discussed above are replaced by more general equa¬ 
tions. The Schrodinger equation is replaced by the well- 
known master equation [17], while the Heisenberg equa¬ 
tion of motion for operators turns into the quantum 
Langevin equation [14,15], which typically takes the form 

d t A=- l -[A,H}-^{[B,A],F- 1 d t B}, (8) 

where {.,.} denotes the anti-commutator, B is the system 
operator that is involved in the coupling to the reservoir, 
F is a noise operator, and 7 a decoherence rate. 

In phase space, the evolution of the Wigner function 
for open quantum systems leads generally to a strictly 
positive Wigner function [18, 19], thus introducing clas- 
sicality [20]. Its dynamical equation takes the form of a 
Fokker-Planck equation [21, 22]. 

To our knowledge, the extension of the Moyal equation, 
i.e., the last of the four dynamical equations, has only 
been addressed in general terms by Ozorio de Almeida 
[23]. In this paper we aim to shed more light on this 
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case by studying the Moyal equation for a Schrodinger 
particle coupled to a thermal reservoir of harmonic oscil¬ 
lators. Our main result is that for Markovian reservoirs, 
in which correlations decay on a short time scale, the 
Moyal equation takes the form 

8 A 

d t A(q,p, t) = {A, H} m - + F(t)d p A , (9) 

where F(t) is a fluctuating force introduced by the reser¬ 
voir, and SA/SF(t) denotes the functional derivative with 
respect to this force. The decoherence rate 7 determines 
the time scale 7 -1 on which information about reservoir 
correlations is lost. 

II. OPEN MOYAL EQUATION 

To derive an example for an open Moyal equation, we 
consider the model proposed by Ford, Kac and Mazur 
[24, 25], in which a single, one-dimensional Schrodinger 
particle is coupled to a reservoir of N harmonic oscilla¬ 
tors. The Weyl symbol of the Hamiltonian is given by 
[32] 

H(q,p,q n ,p n ) = H s + H int +Hr (10) 

Hs =-^-P 2 + V(q) (11) 

2m 

N / 1 x \ 

Hr + Hint — y ] f 2 m Pn ”1” — 9 ) ) > (12) 

n —1 ' n ' 

with k n = m n uin- Here q,p are position and momen¬ 
tum of the system particle, and q n ,p n are the respective 
quantities for the nth reservoir oscillator. The symbol of 
Hr of the reservoir Hamiltonian alone can be found by 
setting q = 0 in Eq. (12). Using the two different repre¬ 
sentations (5), ( 6 ) of the star product, the Moyal bracket 
can be written as 

{A,H} m = ^{H(K) - H(C))A(q,p,q n ,p n ,t). (13) 

It therefore corresponds to linear differential operator on 
phase space that is defined through Eq. (13). It is shown 
in App. A that by a transformation 

A{t) = e t{ '’ ffR}M e 9 if(t) A(t), (14) 

the Moyal equation can be cast into the form 

d t A = {A,H s } M + -K(t)A + F(t)d p A. (15) 
m 

Here F(t) is a function of the reservoir’s phase space 
variables q n ,p n and K(t) is a differential operator acting 
on functions of the reservoir variables. Their definition is 
given in Eqs. (A15) and (A16), but for practical purposes 
only the relation 

K(t)F(t') = —C(t — t') (16) 

N 

C{t ) = y k n cos(o j n t) 

n= 1 


is relevant. The term proportional to F(t) in Eq. (15) 
could be introduced by modifying the system’s potential 
energy as 

V(q)->V(q)-qF(t). (18) 

Therefore, F(t) can be interpreted as a time-dependent 
homogeneous force acting on the system. Because it de¬ 
pends on the reservoir variables, it will fluctuate with the 
state of the reservoir. We therefore can interpret F(t) as 
the phase-space symbol of a fluctuating random force. It 
is the equivalent of noise operators that appear in the 
quantum Langevin equation (8) and therefore may be 
called a noise (Weyl) symbol. 

As K(t) corresponds to a derivative operator, it does 
not represent a phase space symbol, but rather a phase- 
space super-operator that maps the symbol of an oper¬ 
ator to another symbol. We will see below that K(t ) 
describes the dissipation in the system associated with 
the fluctuating random force. 

The function C(t) of Eq. (17) can be interpreted as a 
force correlation function. It is shown in App. B that for 
a thermal reservoir of temperature T hu> n /kB one has 

(F{t)F(t'))=k B TC{t-t'). (19) 


III. MARKOVIAN MOYAL EQUATION 

In many cases of interest, the correlation function (17) 
decays on a time scale r that is much shorter than the 
time scales relevant for the evolution of the system. A 
physical example would be an atom as a system that is 
coupled to the quantized electromagnetic field as reser¬ 
voir. This coupling is responsible for spontaneous decay 
of an excited atom, which for optical transitions happens 
on a time scale of nano seconds. The electromagnetic 
correlation function at these frequencies decays on the 
scale of femto seconds. 

When the reservoir correlation function decays quickly, 
it is possible to make the Markov approximation. To 
analyze the Markov approximation in phase space, we 
first rewrite Eq. (15) as an integral equation, 

A{t) = e t{ > ffs}M A(0) + [ dt'e {t - t,){ -' Hs }m 
J 0 

x (^K{t') + F(t!) d p ) A(t'). (20) 

We note that A(t, q,p , q n ,Pn) is a function of both system 
and reservoir variables, while for a system operator the 
initial Weyl symbol A{0,q,p) does not depend on the 
reservoir. By expanding Eq. (20) into a Dyson series one 
can see that A(t) depends on F(t') through integrals of 
the form 

1= f dt'g(t')F(t'), 

Jo 


(17) 


( 21 ) 
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with g(t') being known functions, and convolutions of 
such integrals. The action of K(t) on such expressions 
amounts to 

K(t)I = - [ dt’g(t')C(t - t'). (22) 

Jo 


so that fa (t) = fa ( 0 ) and 

fa(t) = f dt' fa{t')F{t'). 
Jo 

Inserting this into Eq. (30) yields 


( 33 ) 


To perform the Markov approximation, we define the de¬ 
cay rate 


d t /31 = —fa - 7/^1. 
m 


(34) 


1 f°° 

7 = — / dtC(t). 
m Jo 

Then, for times r, we find 


(23) 


which is solved by 

ACf) = /M 0 )e-^ 


fe(0) 

777,7 


(1 ~ e -7 *). 


(35) 


K(t)I 


- f dt"g(t - t")C(t") 
Jo 

-g(t) f dU'C{f) 

Jo 


—m'yg(t) 


—2my 


SI 


(24) 

(25) 

(26) 
(27) 


With this approximation [33], the open Moyal equation 
takes the form (9). In the following section we will verify 
this approximation at the example of a free particle. 


IV. FREE PARTICLE COUPLED TO A 
RESERVOIR 

To illustrate the general framework of the open Moyal 
equation we consider the situation in which the system 
particle is not subject to an external potential, V = 0. In 
this case the Markovian Moyal equation takes the simple 
form 

T) - ft A 

d t A = —d q A - 2 1P —— + F(t)d p A. (28) 
TO SI ( t ) 


A. Canonical Weyl symbols 

We first solve Eq. (28) for the canonical variables, 
where A( 0 ) = q or A( 0 ) = p, by making the ansatz 

A{q,p, t) = /3i (t)p + fa(t)q + fa(t), (29) 

where fa , fa are functions of time only while fa ( t ) = 
(3o(t,q n ,p n ) may also depend on the phase-space vari¬ 
ables of the reservoir. We use bold greek letters to indi¬ 
cate such a dependence. Inserting this into Eq. (28) and 
sorting the result with respect to q and p yields 

=> - 2 ->im (30) 

d t fa = 0 (31) 

d t fa = F(t)fa, (32) 


For the symbols of position (momentum) we have fa (0) = 
0 and /? 2 ( 0 ) = 1 (/3i(0) = 1 and fafa) = 0 ), respectively, 
so that 

A q (t) = q+ -£-(l - e -7 *) + — f dt' (1 - e-^)F(t') 

TO7 TO7 J 0 

(36) 

Ap(t) = e-^p + [ dt' e" 7 *' F{t'). (37) 

Jo 

Here we have adopted the notation that an index at the 
symbol of an operator refers to its initial value, e.g., 
A q (0) = q. It remains to perform transformation (14). 
By using e LRt F{t') = F(t' — t) and approximation (27) 
we obtain 

A q (t) = qe~^ + ^( l-e" 7 *) 

777-7 

1 / 

+- / dt' (1 - e -7 * )F(t' - t) (38) 

m 7 Jo 

Ap(t) = e -7t (p — nri'yq) + I dt' e -7 * F(t' — t ). (39) 

Jo 

We have verified that this solution agrees with the 
corresponding operator-valued solution of the quantum 
Langevin equation. Furthermore, it is shown in App. C 
that the Markovian solution presented here agrees with 
an exact treatment based on a specific correlation func¬ 
tion C(t). 

To turn solutions (38) and (39) into symbols for the 
system particle alone, we have to take the average with 
respect to the reservoir. It is shown in App. B that for a 
thermal reservoir the mean noise force vanishes, ( F(t )) = 
0. Therefore, the open Weyl symbols of position and 
momentum take the form (38) and (39) with F(t) set to 
zero. 

The physical interpretation of this motion is as follows. 
The coupling to many oscillators with different frequen¬ 
cies results in a dissipation of energy from the particle 
into the reservoir. Therefore, the momentum of the par¬ 
ticle is damped on a time scale 7 -1 . The damping of the 
initial position q arises because in the model by Ford, 
Kac and Mazur all oscillators pull the system particle to¬ 
ward the common equilibrium point at q = 0. However, 
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the particle continues to move in its initial direction for 
a time of the order of 7 -1 , which explains the term pro¬ 
portional to 7 -/( 7717 ) in A q (t). 

The reader may have noticed that A p ( 0) = p — myg 
does apparently not fulfill the correct initial conditions. 
However, this is merely a consequence of the fact that the 
Markovian approximation is only valid for times i > r. 
The non-Markovian derivation in App. C shows that the 
term proportional to 7717(7 builds up on the short time 
scale r and then decays on the long time scale 7 ~ 1 . 


B. Variance Weyl symbols 

To fully appreciate the influence of the reservoir one 
also needs to study the uncertainty of the canonical vari¬ 
ables. We therefore consider symbols with initial condi¬ 
tion A( 0 ) = q 2 or A( 0 ) = p 2 by making the ansatz 

A = /3\p 2 + P 2 pq + ^3Q 2 + PiP + A Q + @6- (40) 

As before, we sort the terms with respect to the power of 
q and p. For four of the six coefficients this can be done 
exactly as for the canonical symbols, leading to 0 3 (t) = 
/3 3 (0) and 


jSa(*) = l-e-T*) 


77Z7 

nt 


The remaining two coefficients obey 


1 


dti 3 i = — #2 - 27 


- 2^1 


SF(t)J 


1 


d t p4 = + - /3 5 - 2 7 


(41) 


(3 5 (t)= [ dt'p 2 (t')F(t') (42) 

Jo 

(3 6 (t)= [ dt' ^(t')F(t'). (43) 

Jo 


-2 


SF(t)J 


(44) 

(45) 


Using Eq. (43) leads to 

Pi{t) = 2/ dt'e'V-Vp^Fit') 

Jo 

+ — r dt' f dt" e 7 ^- t )/3 2 (t ,, )F(i"). (46) 
m J 0 Jo 

The functional derivative of (3^ can then be simplified by 
noting that, for an arbitrary function g{t,t"), 

f dt' f dt"g(t',t")= f dt" f dt'g(t',t"), (47) 
Jo Jo Jo Jt" 

to obtain 


dtfh = —P 2 ^ 27/3i, 
m 


(48) 


which is solved by 

Pi{t) = e-^/MO) + ^(1 - e -7 *) 2 . (49) 

m z 7 z 

Applying transformation (14) to Eq. (40) and using the 
expressions for /3* (i = 4, 5, 6) yields 

A(t) = A{t) - m^q(2{p - m^q)j3i + qf} 2 + At), (50) 

where F(t') is replaced by F(t — t') everywhere. 

Because ( F(t )) = 0, the averaged form of Eq. (50) 
takes the form 

A(t) = Pxp 2 + q 2 (/3 3 - 7717^2 + 2m 2 7 2 /3i) 

+ pq(fe -2m7A) + </3e(0>- ( 51 ) 

The mean value of (3q can be evaluated using Eqs. (43), 
(46), and (19). Within the Markovian approximation, 
the result is given by 

(/3 6 (f)} = 2k B Tm r ) f dt'/3i(t') (52) 

Jo 

= mk B T{p 1 ( 0)(l-e- 27 ‘) 

+ & ( 0 )( 2 7 t -e-» + 4^‘-3)i. 

7 z m z 1 

It remains to apply the initial conditions, /3i(0) = 1 for 
A p 2 (t) and 0 3 (0) = 1 for A q 2 (t), with all other coefficients 
being zero. Putting everything together we arrive at the 
open Weyl symbols 


A q 2 {q,p,t) = A 2 + q 2 (l-e 7t ) 2 

+ (2 7 t - 3 - e~ 2lt + 4e -7 *) (54) 

777-7 z 

A p i{q,p,t) = A 2 + mk B T (l - e _27t ) + e~ 2lt m 2 ^ 2 q 2 . 

(55) 

We have confirmed that this solution agrees with the 
operator-valued solution derived from the quantum 
Langevin equation. The physical interpretation is as fol¬ 
lows. The term 2k B Tt/ (7717) in A q 2 corresponds to a 
diffusion of the particle position. For solutions of the 
diffusion equation d t f(t,q) = Dd q f(t,q), the variance of 
position increases as 2Dt for sufficiently large times. This 
implies that the coupling to the reservoir can be linked 
to a diffusion coefficient D = k B T/ (7717). The Einstein- 
Smoluchowski relation D = gk B T then implies that the 
mobility of the system particle is given my g = 1/(7717). 

The term proportional to mk B T in A p 2 describes the 
thermalization of the system particle through its cou¬ 
pling to the reservoir. Keeping in mind that (Ekm) = 
(A P 2 )/ (2m) one can see that for times f > 7 -1 the kinetic 
energy of the reservoir approaches k B T / 2 , confirming the 
equipartition theorem. Finally, the terms proportional to 
q 2 are a consequence of the dragging towards the origin 
in the Ford-Kac-Mazur model that we discussed above. 
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V. DISCUSSION AND CONCLUSION 

In the previous sections we derived the open Moyal 
equation (9), which describes the dynamics of open quan¬ 
tum systems coupled to a Markovian reservoir, and illus¬ 
trated its features at the example of a free particle. The 
significance of the open Moyal equation is not so much in 
the specifics given here, but in the possibility to extend 
this model to describe open systems of larger interest. 

This is the same situation as with master equation and 
quantum Langevin equation ( 8 ). Both equations have 
been derived from specific models, but have been general¬ 
ized to describe the influence of various reservoir-induced 
effects on a large variety of systems. Examples include 
spontaneous emission, thermal excitation, spin dephas¬ 
ing, vibrational relaxation in molecules, and lossy optical 
cavities [14, 15]. 

In practice, most researchers do not derive reservoir 
properties from first principles but rather pick an ad hoc 
model to include the effect of a reservoir. This choice is 
not arbitrary but has to obey general principles. Posi¬ 
tivity and trace preservation of the density matrix limit 
Markovian master equations to the celebrated Lindblad 
form [17]. Likewise, noise operators F(t) and decoher¬ 
ence rate 7 in quantum Langevin equations are related 
through a fluctuation dissipation theorem. 

The same is true for the open Moyal equation. A rather 
trivial extension of the free particle discussed above is a 
particle under the influence of a constant external force 
Fq. This can be accomplished by including a linear po¬ 
tential V(q) = —qF 0 in Eq. (A2). Using the techniques 
of Sec. IV it is not hard to see that the only change in 
solutions (38), (39) is then to replace F(t) by F(t) + Fq. 
The Weyl symbols of position and momentum are then 
modified according to 


A q {t) -a A q (t) + ^ (t - j-\l - e "*)) (56) 

777-7 

A p (t) —> A p (t) H- (1 — e lt ). (57) 

7 

For large time t this describes a particle moving with 
drift velocity Fq /(my), which confirms the result for the 
mobility p found in Sec. IV B. Hence, the system’s re¬ 
sponse to a linear perturbation is linked to the diffusion 
of q through the Einstein-Smoluchowski relation, which 
is one of the earliest examples of a fluctuation dissipation 
theorem. 

The model presented here can be readily extended in 
several ways. An obvious extension is to combine the 
model by Ford, Kac and Mazur with more complicated 
external potentials V(q). An interesting question arises 
for nonlinear potentials of Kerr type, where the exact 
solution for Weyl symbols exhibits a singularity [13]. Be¬ 
cause dissipation often has a moderating effect, it may 
be possible that coupling to a reservoir may eliminate 
this singularity. One may also consider the effect of a 
reservoir on a system that consists of two or more inter¬ 


acting particles, for instance for two photons interacting 
via cross Kerr modulation [26]. 

Another extension would be to generalize the quadratic 
coupling between the position of system and reservoir 
particles in the Ford-Kac-Mazur model. One way would 
be to introduce a distribution of equilibrium positions for 
the reservoir oscillators and to average over this distri¬ 
bution. This may eliminate the drag towards the origin 
that we discussed in Sec. IV. One may also introduce a 
different coupling like P^2 n Pn in Eq. (12), which could 
be realized by a series of quantum LC circuits coupled 
through their mutual inductance [27]. In the context 
of quantum Langevin equations, it is possible to extend 
the interaction with a reservoir to nonlinear couplings 
[26, 28, 29], although this is considerably more involved. 
It is conceivable that the same could be accomplished for 
the open Moyal equation. 

Finally one could follow the common practice in the 
field of quantum Langevin equations and simply pick an 
ad hoc model to obtain another open Moyal equation. 
This could be done by replacing p in Eq. (9) by another 
system observable O , and F(t) by random fluctuations of 
the corresponding time derivative d t O. Given that phase 
space methods are an excellent tool to compare classical 
and quantum dynamics, and given that coupling a quan¬ 
tum system to a reservoir generally suppresses quantum 
correlations, the open Moyal equation may therefore pro¬ 
vide a promising method to study how classical behaviour 
emerges in open quantum systems. 
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Appendix A: Derivation of the non-Markovian open 
Moyal equation 

To avoid a cluttered notation, we will define linear op¬ 
erators Li through 


LiA := { A , Hi} M , i = S,R, int. 


(Al) 


The explicit form of the Moyal bracket for the various 
parts of the Hamiltonian is then given by 


Ls= ^d q +^[V(q-i^d p )-V (q + i^d p )} 


N 


Lr + Lint — 'y ] 


—dq n - k„(q n - q)(dp n - d p ) 
TYl n 


(A2) 


(A3) 
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To solve the reservoir part of Eq. (13), it is useful to 
switch to complex coordinates 


1 ( Qn .In \ 

“"“cnW’-V”) 

(A4) 

L-J h , 

(A5) 

V m n uj n 

which yields 


N N , , 


Lint - qdp E kn + iq E /r )h (&*„ 

n — 1 n .— 1 

-5<) 


+ ( A6 ) 

n= 1 V 
N 

Lr — ZC<J n (ck n ^Q;* Q-ndoLn)' (AT) 

n—1 

We now go into an interaction picture with respect to the 
reservoir by setting 

A(q,p,a n: a^) = e tLR A{q 1 p, a n , a* ). (A8) 

The Moyal equation then takes the form 

d t A = e- tL «(L s + L int )e tL «A (A9) 

N N b i 

= Ls A- q d p J2knA+J2^ 

n—1 n—1 *^ 

X ~ d <X-n) etLR 

+ e~ tLR (a n + a* n )e tLR d p )A. (A10) 

Operators of the form 0(t) = e~ tLR O{0)e tLR obey the 
differential equation d t O = —[Lr,0\. It is not hard to 
see that for 0(0) = a„ and 0(0) = d an , this equation 
has the solution 


e~ tLR a n e tLR = e iuJnt a n (All) 

e~ tLR d an e tLR = e~ iUnt d ari . (A12) 

Hence the Moyal equation can be rewritten as 

/ N \ ~ N k 1 

d t A= lLs-qd p ^Tk n \A+Y^ 

\ n— 1 ) n—1 V ^ 

x {jA {e~ iuint d an - e iUnt d a *) 

+ (e iUnt a n + e-“”*a;)5 p }l (A13) 

= {L s - C(0)qd p + Kq + F(t) d p }A, (A14) 

with 

Hi) = E + e-^a* n } (A15) 

n—1 ^ ^ 

N lc 1 

K(t ) = - E {e~ iuint d an + } (A16) 

„_i \AnWn 


To perform the Markovian approximation it is useful to 
reformulate Eq. (A14) by setting A(t) = exp (qK(t))A(t). 
Because [K(t),K(t')] = 0 we can use the Baker- 

Campbell-Hausdorff formula and Eq. (16) to transform 
the open Moyal equation into 

d t A = {e~ qK ^L s e qK ^ - C(0)qd p 


+ e - qKW F(t)e qK{t) d p } A 

(A17) 

{L s + —K{t) - C(0)qd p 
m 


+ (F(t)-q[K(t),F(t)])d p }A 

(A18) 

{Ls H- E(t) A F(t) d p } A, 

(A19) 


which is Eq. (15). 


Appendix B: Averaging over the reservoir degrees of 
freedom 


To obtain Weyl symbols that only depend on the sys¬ 
tem variables, one has to evaluate their mean value with 
respect to the reservoir degrees of freedom. Generally, 
the mean value of a symbol A(q,p) in phase space is ex¬ 
pressed through 

(A) = J dqdpA(q,p) w (q,p). (Bl) 


For a harmonic oscillator in thermal equilibrium, the 
Wigner function takes the form [30] 

W = W ^re- H(q ' p)/{E) (B2) 

— u -M«l 2 /(e) (-po\ 

“2 n(E) e 


with (E) = coth(Huj/(2kBT)) the mean en¬ 

ergy and H(q,p) = p 2 /(2m) + \muj 2 q 2 the Weyl 
symbol of the Hamiltonian. This result can be 
derived using the thermal density matrix p = 
exp[— i?/(fcsT)]/Trexp[— H/{ksT)] for the harmonic os¬ 
cillator and Mehler’s formula (see Sec. 10.13 of Ref. [31]). 
If we assume that the reservoir oscillators are thermal- 
ized, one can easily evaluate that (|a n | 2 ) = (E n ) / (fruj n ) 
and (q„) = (ck 2 ) = 0. Eq. (A15) then implies ( F(t )) = 0 
and 


<F(i)F(t'))=E^n 

n 


(En) 

tuxJn 


COS (w n {t — t')). 


(B4) 


For a hot reservoir, where ( E n ) Huj n , one has (E n ) « 
ksT, from which Eq. (19) follows. 


Appendix C: Solution for canonical observables 
without Markov approximation 

In this section we repeat the calculations of Sec. IV A 
without making the Markov approximation (27). This 
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amounts to solving Eq. (28) with — 2 q 6A/8F(t) replaced 
by K(t)A/m. By making the same ansatz as in Sec. IV A 
we obtain the same equations and solutions for 02 and 
0 3 , but is determined by 


which fulfills J^° C(t)dt = my. The parameter 7 should 
correspond to the decay rate associated with dissipation, 
while t = 1 /T corresponds to the short time scale on 
which C(t) decays. 


"101 = 02 + AT0 3 . 

Using Eq. (16) this can be transformed into 


(Cl) 


TO01=02(0)-/ dt'0i(0)<C(i-f'). (C2) 


Jo 

This equation can be solved using a Laplace transforma- 
tion g(t) -> g(s), 


If Si denotes the poles of C(s), then the inverse 
Laplace transformation contains time-dependent expo¬ 
nentials e tSi multiplied by constant factors. For T> 7 , 
we can make a separate Taylor expansion of the poles 
in the exponent and the constant terms around T = 00 . 
The poles Sj are then approximately given by —7 and 
—T. After applying transformation (14) we obtain, to 
leading order in 1/T, A q (t) of Eq. (38) as well as 


0iO) = 


S + m C ( S ) 


0i(O) + —02(0) . (C3) 


As an example for a reservoir memory function we con¬ 
sider [34] 


C(t) = m-/Te- rt 
myT 


C(s) = 


s + r’ 


(C4) 

(C5) 


A p = e lt p — m r yq(e 11 — e rt ) + 



dt' e-i*'F{t-t'). 

(C 6 ) 


For times t^> r this solution agrees with the Markovian 
result (39). It only differs by the short-living term pro¬ 
portional to e _rt , which ensures that the correct initial 
conditions are fulfilled. 
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The Hamiltonian operator takes the same form, with the 
symbols q,p,q n ,p n replaced by the respective operators 
because S[<?(<))] s= g(q) for any operator g(q) that is a 
function of q only. A similar statement holds for operators 
of the form g(p). 

While 8I/8F(t') = g(t') for all points 0 < t! < t, one 
has to be careful at the boundary of the integral I. The 
factor of 2 in Eq. (27) is a consequence of the follow¬ 
ing procedure. We rewrite I by replacing the integration 
boundary t by 00 and the function g(t') by g(t')9(t — t'), 
with 9 the step function. Then SI/SF(t) = g(t)8( 0). Set¬ 
ting 0(0) = 1/2 yields Eq. (27). 

Eq. (17) suggests that dtC(0) = 0. Although example 
(C4) does not obey this relation, we have verified that 
in the limit F > 7 it yields the same result as more 
complicated examples that do fulfill it. 
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